Topological Catastrophe and Isostructural Phase Transition in Calcium 
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We predict a quantum phase transition in fee Ca under hydrostatic pressure. Using density 
functional theory, we find at pressures below 80 kbar, the topology of the electron charge density is 
characterized by nearest neighbor atoms connected through bifurcated bond paths and deep minima 
in the octahedral holes. At pressures above 80 kbar, the atoms bond through non-nuclear maxima 
that form in the octahedral holes. This topological change in the charge density softens the 
elastic modulus of fee Ca, while C44 remains unchanged. We propose an order parameter based on 
applying Morse theory to the charge density, and we show that near the critical point it follows the 
expected mean- field scaling law with reduced pressure. 
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The theory and characterization of solid-solid isostruc- 
tural phase transitions has been an area of intense 
experimental and theoretical research since they were 
described in 1949 [T. What sets these transitions 
apart from others is that while crystalline proper- 
ties change-sometimes discontinuously-crystallographic 
structure does not, making these transitions purely elec- 
tronic in nature. Both first order and continuous tran- 
sitions are known to exist. A well known example of 
the former is the 7 ^ a transition in fee Ce [1 , while a 
magneto- volume instability in fee Fe is an example of the 
later |2j. 

In the case of Fe, the traditional local order parame- 
ter, //, gave conflicting results computationally, with the 
transition appearing as first [3 or second [2 order. It was 
not until the transition was interpreted as a topological 
change of the elctronic charge density, p(r), that the un- 
derlying nature of the phase change became apparent [4^ . 
While this approach is promising, no general correlations 
have been drawn between isostructural phase transitions 
and p(r). We use this approach to predict an isostruc- 
tural phase transition in calcium under hydrostatic pres- 
sure. We then apply Morse theory to the charge density 
and propose an order parameter for the transition. We 
show that this order parameter scales in the appropri- 
ate way with reduced pressure near the quantum critical 
point. 

One of the advantages of casting theories of isostruc- 
tural phase transition in terms of p{f) is that it is a quan- 
tum observable. Though p(r) is often calculated by way 
of first principles methods (e.g. density functional theory 
(DFT)), the total charge density can also be measured 
via X-ray diffraction techniques, and the spin-polarized 
charge density can be determined using spin-polarized 
neutron diffraction. As it is also known that all ground 
state properties depend on the charge density [5 , it seems 
appropriate to seek relationships between the structure 



* Electronic address: |trjones@mines.edu| 



of p(r), changes to that structure, and corresponding 
changes to properties. We note that previous work has 
proposed that DFT could be used as a general framework 
for analyzing quantum phase transitions in electronic sys- 
tems [6]. 

The electron charge density is a 3 dimensional scalar 
field. Morse theory tells us the charge density can be 
characterized by its rank 3 critical points (CPs), the 
points at which the gradient of the field vanishes. There 
are four kinds of rank 3 CP in a three-dimensional space: 
a local minimum, a local maximum and two kinds of 
saddle point. These are conventionally denoted by an in- 
dex, which is the number of positive curvatures minus the 
number of negative curvatures. For example, a minimum 
CP has positive curvature in three orthogonal directions; 
therefore it is called a (3, +3) CP. The first number is 
the dimensionality of the space, and the second number 
is the net number of positive curvatures. A maximum is 
denoted by (3, -3), since all three curvatures are negative. 
A saddle point with two of the three curvatures negative 
is denoted (3, -1), while the other saddle point is a (3, 
+1) CP. 

Through extensive studies of molecules and crystals, 
Bader and Zou showed that it was possible to correlate 
topological properties of p{r) with elements of molecular 
structure and bonding |7|E]- In particular, a bond path 
was shown to correlate with the ridge of maximum charge 
density connecting two nuclei. The existence of such a 
ridge is guaranteed by the presence of a (3, -1) CP. As 
such, this CP is often referred to as a bond CP. Other 
types of CPs have been correlated with other features 
of molecular structure. A (3, +1) CP is topologically 
required at the center of ring structures. Accordingly, it 
is designated a ring CP. Cage structures must enclose a 
single (3, +3) CP and are given the name cage CPs. The 
locations of the atomic nuclei always coincide with a (3, 
-3) CP. Hence, it is often called an atom or nuclear CP. 

In this study, the charge density of fee Ca was calcu- 
lated at various lattice constants and strains using the 
Vienna ab-initio Simulation Package (VASP) version 4.6 
[9]. All calculations were performed with the Perdew- 
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Burke-Ernzenhof generalized gradient corrections [10], 
and while we did investigate the effect of unpaired spins, 
Ca was found to be non- magnetic. As such, the results 
reported here are for spin-restricted PBE calculations. 
Forces for the calculation of phonons were found using 
5x5x5 supercells in VASP. Frozen phonon calculations 
were performed using fropho [11 . 

Calcium at zero pressure and temperatures below 721 
K crystallizes into an fee phase. At 300 K it remains fee 
up to a pressure of 195 kbar, at which point it trans- 
forms to bcc [12 . In this study, we focus attention on a 
region around the low pressure fee phase and find that it 
contains an isostructural transition. 

Most nonmagnetic, monoatomic fee metals share the 
ground state charge density topology of the prototype, 
Cu [T?, with each atom connected to its 12 nearest neigh- 
bors by single bond paths. The ground state charge den- 
sity of Ca, however, is characterized by bifurcated bond 
paths, where bond point pairs are separated by a (3,+l) 
CP (Fig. [T]). Like the prototype, however, there is a 
(3, +3) CP at the center of the octahedral hole. These 
features can be seen in the plot of p{f) in the (200) plane 
shown in the center of Fig.[T] The position of the plane in 
the unit cell is shown in the left of Fig. [l] Minima in this 
plot are shaded black, while the maxima are unshaded. 
Ca atoms lie on the edges of the plot. 

As the pressure is increased towards the fcc-bcc tran- 
sition pressure of 195 kbar, the bifurcated bond paths 
separate by moving towards the nearest octahedral cage 
CP. The four bond points nearest the octahedral hole 
in the center of the plot move towards that cage CP, 
while the other bond points move in the opposite direc- 
tions, towards the octahedral holes at the corners of the 
plot (shown in the right of Fig. [T]). By symmetry, all of 
the CPs move at the same rate. Thus, the four bond 
points nearest the octahedral cage point in the center of 
Fig. [l] will meet the cage point simultaneously, along with 
the four bond points on the two perpendicular planes, 
(020) and (002). Ah 12 CPs coalesce with the octahe- 
dral cage point when the external pressure reaches 79 
kbar This topological catastrophe produces a non- 
nuclear maximum at the center of the fee unit cell, the 
octahedral hole. 

The maxima and bond paths for this high pressure 
topology are shown in Fig. [2j Here, the larger gray 
spheres represent the Ca atoms, while the smaller white 
ones represent the non-nuclear maxima. Inspection of 
this figure reveals that all Ca atoms are now bound to 
the 12 nearest neighbor and 6 second neighbor atoms 
through the non-nuclear maxima located in the octahe- 
dral holes of the fee lattice. The contour plot in the right 
of Fig. [2] shows the contours of p{r) in the (200) plane. 
This is the same crystallographic plane shown in Fig. [l] 
with the same shading for minima and maxima. Com- 
paring the two contour plots, it is clear that the CP in the 
octahedral hole has changed from a minimum to maxi- 
mum while the CP on the atom-atom line remains a ring 
point. There are no further changes in the topology as 



the pressure is increased to 195 kbar. 

A key result of this work is the discovery of a pressure- 
induced change to the topology of p{f) that is not ac- 
companied by a change in the crystallographic structure 
of fee Ca. This transition occurs at a pressure of 79 
kbar when the lattice constant is c^5.0 A. The transi- 
tion is the result of a topological catastrophe, through 
which the cage point at the center of an octahedral hole 
is transformed to a non-nuclear maximum. 

Because there is an abrupt change to calcium's charge 
density topology, we expect that properties measuring 
the response of p{f) to an external perturbation will also 
change through this transition (e.g. the elastic properties 
of the material). To begin, we investigated the material's 
phonons, and while most stiffen with increasing pressure, 
the transverse shear wave propagating in the [110] direc- 
tion with polarization in the [110] direction softens after 
the topological catastrophe. The stiffness of this mode is 
given by the C single crystal elastic modulus. 

To understand the softening, consider both single crys- 
tal elastic moduli, C44 and C^ These were calculated us- 
ing VASP by employing a finite displacement approxima- 
tion [14]. At the equilibrium lattice constant, these cal- 
culated single crystal elastic constants are in good agree- 
ment with experimental values. We found 0.046 Mbar 
and 0.16 Mbar for C and C44 respectively, compared to 
the experimental values of 0.047 and 0.16 [15]. In the low 
pressure phase, both C44 and C increase with applied 
pressure, while in the high pressure phase C44 increases 
with applied pressure and C decreases (Fig. |3|, with C 
achieving its maximum value at the topological catastro- 
phe. This can be understood in terms of screening of 
the nuclear charges. The strains resisted by C44 and C 
are labeled A and B respectively in Fig. [4| Strain A has 
non- vanishing Cxy^ while strain B has €xx = ~^yy Non- 
nuclear maxima in the high-pressure phase screen the 
restoring force for strain B and soften the elastic con- 
stant C ^ while no such screening is present to soften C44 

We have found that the topological catastrophe of p{f) 
is the result of changes in the occupancy of calcium's one 
electron orbitals. Specifically, calculations on 55 atom 
clusters of fee Ca at various lattice constants using the 
Amsterdam Density Functional (ADF) package fT6| re- 
vealed that at the lattice constant corresponding to the 
catastrophe, filling an orbital comes at the expense of 
filling a T2g orbital. The latter is composed of s charac- 
ter on the second coordination sphere, while the former 
contains dx2-y2 and 6.^2 character on the central atom, 
which contributes density to the octahedral hole of the 
fee lattice. When the E^ is filled, the extra density in the 
octahedral holes due to the anti- nodes on the dx2-y2^ 6.^2 
orbitals causes the topological change of p{f). 

We propose using the distance s between the bond CPs 
and the cage CP as a convenient order parameter for 
the transition. As the pressure is increased towards the 
critical pressure, Pc, the CPs around the octahedral hole 
move towards the center of the hole, as shown in Fig. [l] 



FIG. 1: Atoms and bond paths in the low pressure topology of fee Ca are shown left. Contours of the total charge density 
are shown for the (200) plane (center). The shaded regions represent minima in p(f). The bifurcated bond paths and arrows 
indicating their motion under pressure are shown far right. 




FIG. 2: The atoms and bond paths in the high pressure 
topology of fee Ca (left). Contours of the total charge density 
are shown for the (200) plane (right). The shaded regions 
represent local minima in p{r). 




FIG. 4: 
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Strains with only non-vanishing (A) txy and (B) 
yy. The dashed line indicates the strained square. 



Non-nuclear maxima in the high-pressure phase are denoted 
by open circles, while filled circles denote Ca atoms. 
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FIG. 3: Plot of C44 (crosses) and C^ (circles) in Mbar as a 
function of pressure in kbar. The C^ points are also connected 
by a dashed line to help guide the eye, and the lattice constant 
corresponding to the topological change is indicated with a 
vertical dashed line. 



until they coalesce at the quantum critical point (QCP) 
with P=Pc where s vanishes. On the basis of numerical 
studies near the QCP, we show that s follows the mean- 
field scaling law with reduced pressure p= \P — Pd/Pc 
Using the well-known mapping between quantum sys- 
tems at zero temperature in d-dimensions and classical 



systems in d -\- 1-dimensions, we anticipate that in 3d a 
scalar order parameter near the QCP should scale with 
reduced pressure diS s ^ ps with the mean-field criti- 
cal exponent S = 3 taken. Fig. [5] shows DFT results of 
\np/\ns near the QCP. 

A very fine charge density mesh was required to 
achieve the desired accuracy. However, our results re- 
mained unchanged at grid densities beyond 200 points 
per Angstrom. As such, the results reported here were 
calculated using 200 points per Angstrom mesh. We find 
that 1/S reaches a value of 0.35 dit p 8e-5, and converges 
to 0.33 upon linear extrapolation to p = 0. (Quadratic 
and cubic extrapolations were also performed using the 
6 points nearest p = 0, leading to 0.34 and 0.33, respec- 
tively.) 

In summary, we have used the topology of p{f) to pre- 
dict an isostructural phase transition in fee calcium. We 
found that fee calcium adopts different topologies of p{f) 
when the applied hydrostatic pressure is above and be- 
low 79 kbar. We predict that this topological change in 
the charge density softens the C elastic modulus of fee 
Ca, while C44 remains unchanged. We propose an order 
parameter based on applying Morse theory to the charge 
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density, and we show that near the critical point it fol- 

6 . • ' " lows the appropriate mean-field scaling law with reduced 
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FIG. 5: Plot of In 5/ In p vs. p. Extrapolation to p = gives 
a limiting value of 1/^ = 0.33, consistent with the mean-field 
scaling law. 
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